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Primordial scalar perturbations may be generated when complex conformal scalar field 
rolls down its negative quartic potential. We begin with the discussion of peculiar infrared 
properties of this scenario. We then consider the statistical anisotropy inherent in the model. 
Finally, we discuss the non-Gaussianity of scalar perturbations. Because of symmetries, the 
bispectrum vanishes identically. We present a general expression for the trispectrum and 
give its explicit form in the folded limit. 
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§1. Introduction and summary 

Primordial scalar perturbations in the Universe are approximately Gaussian and 
have approximately flat power spectrum.^) The first property suggests that these 
perturbations originate from amplified vacuum fluctuations of weakly coupled quan- 
tum field(s). The flatness of the power spectrum may be due to some symmetry. 
The best known is the symmetry of the de Sitter space-time under spatial dilatations 
supplemented by time translations. This is the approximate symmetry of the inflat- 
ing Universe, which ensures approximate flatness of the scalar spectrum generated 
by the inflationary mechanism.'^) Inflation is not the only option in this regard, 
however. Indeed, the flat scalar spectrum is generated also in the scalar theory with 
negative exponential scalar potential in flat space-time^-* (see also ref.^-*). Its equa- 
tion of motion is invariant under space-time dilatations supplemented by the shifts 
of the field. This symmetry remains approximately valid in slowly evolving, e.g., 
ekpyrotic^-* or "starting" Universe, hence the flatness of the resulting perturbation 
spectrum in these models. It is worth noting that there are other mechanisms ca- 
pable of producing flat or almost flat scalar spectrum.^-*'^^ In some cases, there is 
no obvious symmetry that guarantees the flatness, i.e., the scalar spectrum is flat 
accidentally. 

In search for alternative symmetries behind the flatness of the spectrum, one 
naturally comes to conformal invariance.^^^' In the scenario of ref.,^^^ it is supple- 
mented by a global symmetry. The simplest model of this sort has global symmetry 
C/(l) and involves complex scalar field 0, which is conformally coupled to gravity 
and for long enough time evolves in negative quartic potential 

V{ct>) = -/i2|c/>|4 . (M) 

The theory is weakly coupled at /i < 1. One assumes that the background space-time 
is homogeneous, isotropic and spatially flat, ds'^ = a'^{r]){dri'^ — dx^). Then, due to 
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conformal invariance, the dynamics of the field xi'n^^) — '^4' is independent of the 
evolution of the scale factor and proceeds in the same way as in Minkowski space- 
time. One begins with the homogeneous background field Xc(??) that rolls down 
the negative quartic potential. Its late-time behavior is completely determined by 
conformal invariance, 

where r]^ is an arbitrary real parameter, and we consider real solution, without loss 
of generality. As we review in section 2, at early times the linear perturbations about 
this background oscillate in conformal time as modes of free massless scalar field, 
while at late times the perturbations of the phase 

9 = V2Avg (j) 

freeze out. Somewhat unconventional normalization of the phase is introduced for 
future convenience. At the linear level, their power spectrum is flat. 



h 

2^ 



(1-3) 



As discussed in ref.,^"^ this property is a consequence of conformal and global sym- 
metries. 

The scenario proceeds with the as- 
sumption that the scalar potential V{(j)) 
has, in fact, a minimum at some large 
value of 1 1^1, and that the modulus of 
the field (f) eventually gets relaxed to the 
minimum, see figure 1. The simplest 
option concerning further evolution of 
the perturbations 59 is that they are 
superhorizon in the conventional sense 
by the time the conformal rolling stage 
ends. We proceed under this assump- 
tion. The phase perturbations remain 
frozen out,*) and their power spectrum 
remains flat. At some much later cos- 
mological epoch, the perturbations of 
the phase are converted into the adia- 
batic scalar perturbations by, e.g., the 
curvaton mechanism^^^ (in that case 9 
is a pseudo-Nambu-Goldstone curvaton, 

and reprocessing occurs as discussed in ref.^^-*) or modulated decay mechanism. ^^)' 
In either case, the power spectrum is not distorted, so the resulting adiabatic per- 
turbations have flat primordial power spectrum. If conformal invariance is not exact 

*' For contracting Universe, this property of superhorizon modes holds if the dominating matter 
has stiff equation of state, w > 1. This appears to be necessary for the viabihty of the bounce scenario 
anyway, see the discussion in refs.^^'' 




Fig. 1. The scalar potential. Bullets show the 
evolution of the scalar field. Arrows at the 
end point at the bottom of the potential 
indicate that there are perturbations of the 
phase. 
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at the rohing stage, the scalar power spectrum has small tilt, which depends on both 
the strength of the violation of conformal invariance and the evolution of the scale 



A peculiar property of the model is that the modulus of the rolling field also 
acquires perturbations. At late times, modes of the modulus (i.e., radial direction) 
have red power spectrum (see section 2.2 for details). 



One consequence is that there are perturbations of the energy density with red 
spectrum right after the conformal rolling stage, but before the modulus freezes 
out at the minimum of V{(f)). These are not dangerous, provided that the energy 
density of the field (j) is small compared to the total energy density at all times before 
the modulus freezes out, i.e., the cosmological evolution is governed by some other 
matter at that early epoch. In this paper we assume that this is indeed the case. 

The second consequence is that the infrared radial modes interact with the 
perturbations of the phase, and in principle may have strong effect on the latter. 
This is one of the issues we address in this paper. We show that to the linear order 
in h, the infrared effects can be absorbed into field redefinition, so there is no gross 
modification of the results of the linear analysis due to the effect of the infrared 
modes. 

The large wavelength modes of 5\(p\ are not entirely negligible, however. The 
modes whose present wavelengths exceed the present Hubble size H^^ induce sta- 
tistical anisotropy in the perturbations of the phase 69, and hence in the resulting 
adiabatic perturbations: the power spectrum of the adiabatic perturbation ( has the 
following form. 



The first non-trivial term, linear in h, is free of the infrared effects; Wij is a traceless 



vector, k = k/A;, and ci is a constant of order 1 whose actual value is undetermined 
because of the cosmic variance. In the last term, u is some unit vector independent 
of Wij , and the positive parameter C2 is logarithmically enhanced due to the infrared 
effects. This is the first place where the deep infrared modes show up. Clearly, their 
effect is subdominant for small h. 

The statistical anisotropy encoded in the last term in (1-5) is similar to that com- 
monly discussed in inflationary context, and, indeed, generated in some concrete 
inflationary models:^^' it does not decay as momentum increases and has special 
tensorial form (ku)^ with constant u. On the other hand, the first non-trivial term 
in (1-5) has the general tensorial structure and decreases with momentum. The lat- 
ter property is somewhat similar to the situation that occurs in cosmological models 
with the anisotropic expansion before inflation. Overall, the statistical anisotropy 
(1-5) may be quite substantial, since there are no strong bounds on h at least for the 
modulated decay mechanism of conversion of the phase perturbations into adiabatic 
ones. 



factor at the rolling stage. 




(1-4) 




(1-5) 
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The non-linearity of the scalar potential gives rise to the non-Gaussianity of the 
perturbations of the phase 59, and hence the adiabatic perturbations in our scenario, 
over and beyond the non-Gaussianity that may be generated at the time when the 
phase perturbations get reprocessed into the adiabatic perturbations. In view of the 
result outlined above, this non-Gaussianity is not plagued by the infrared effects 
at the first non-trivial order in h. Therefore, the gradient expansion of the effective 
background is useless for the study of the non-Gaussianity, and we have to perform a 
hard-core calculation. Because of the symmetry 9 — t- —0, the bispectrum of the phase 
vanishes, while the general expression for the trispectrum is rather cumbersome. It 
simplifies in the folded limit (according to the nomenclature of ref.^^^); the explicit 
form of the trispectrum in this limit is given by eq. (5-16). 

The paper is organized as follows. In section 2 we review the linear analysis of 
the model. In section 3 we study the effect of infrared modes of the modulus 
on the perturbations of the phase 69 at the leading and subleading orders of the 
gradient expansion, and to the linear order in h. Statistical anisotropy is analysed 
in section 4. Non-Gaussianity is considered in section 5. We conclude in section 6. 

§2. Linear analysis 

At the conformal rolling stage, the dynamics of the scalar field is governed by 
the action 

g^^-^^^*^A+■^ct>*(|>-V{<t>) . 

where the scalar potential V{(p) is negative and has conformally invariant form (1-1). 
In terms of the field x = CL(j), the field equation is 

v'"'d^d,X-2h''\xW = 0. (2-1) 

Spatially homogeneous background approaches the late-time attractor (1-2). 

2.1. Perturbations of phase 

At the linearized level, the perturbations of the phase and modulus of the field 
(j) decouple from each other. Let us begin with the perturbations of the phase, or, 
for real background (1-2), perturbations of the imaginary part X2 = Im x/V^- They 
obey the linearized equation, 

(5x2)" - didi 6x2 - 2h\l <5x2 = , (2-2) 

where prime denotes the derivative with respect to r]. Let k be conformal momentum 
of perturbation. An important assumption of the entire scenario is that the rolling 
stage begins early enough, so that there is time at which 

k{ri^ - r?) > 1 . (2-3) 

Since the momenta k of cosmological significance are as small as the present Hubble 
parameter, this inequality means that the duration of the rolling stage in conformal 
time is longer than the conformal time elapsed from, say, the beginning of the hot 
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Big Bang expansion to the present epoch. This is only possible if the hot Big Bang 
stage was preceded by some other epoch, at which the standard horizon problem is 
solved; the mechanism we discuss in this paper is meant to operate at that epoch. 
We note in passing that the latter property is inherent in most, if not all, mechanisms 
of the generation of cosmological perturbations. 

Equation (2-2) is exactly the same as equation for minimally coupled massless 
scalar field in the de Sitter background. For future reference, we write its solution 
in the following form: 

SX2{^,V) = J (2;r)3/2V2fe {^xi~HK^,v)A + h.c.) . 

Here 

5x^2^ (k, X, T]) = e^kx-ifer,. . p^j^^ n,-ri), (2-4) 

with 



F{k,0 = -^^k^H',/m) (2-5) 
and -^3^2 Hankel function. At early times, the mode oscillates, 

,5x^-^(k,x,r/)=e^'^--^^V (2-6) 

so and A^^ are annihilation and creation operators obeying the standard commu- 
tational relation, [4^,4]^,] = (5(k — k'). As usual, we assume that the field 6x2 is 
initially in its vacuum state. 

At late times, when A;(?7* — 77) ^ 1, the perturbations of the phase are time- 
independent. 

This expression describes Gaussian random field (cf. ref.^^^) whose power spectrum 
is given by (1-3). 

The phase perturbations can be converted into adiabatic ones by at least two 
mechanisms. One operates if is actually pseudo-Nambu-Goldstone field that lands 
at a slope of its potential. This mechanism produces the non-Gaussianity of local 
form in the adiabatic perturbations. For generic values of the phase at landing, 
9c ^ 7r/2, non-observation of the non-Gaussianity^^ implies Vse < 10"*^, cf. ref.,^^) 
so that the correct scalar amplitude is obtained for 

h < 10-2 . (2-8) 

Generally speaking, such a constraint is not characteristic of the alternative, mod- 
ulated decay mechanism. ^^-''^^^ In that case, if the relevant mass or width depends 
linearly on 9, the resuling non-Gaussianity parameter is fairly small (see refs.^^^'^^^ 
for details), Jnl ~ 1, in comfortable agreement with the existing limit. 
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2.2. Perturbations of modulus 

Let us now consider the radial perturbations or, with our convention of real 
background Xc, perturbations of the real part xi = Re x/V^- At the linearized 
level, they obey the following equation at conformal rolling stage, 

{5xi)" - m 6x1 - Qh^ldxi = (dxi)" - didi 6x1 - _ ^y ^Xi = . 

Its solution that tends to properly normalized mode as />:(ry* — ry) — )• oo is 



where Sk, is another set of annihilation and creation operators. At late times, 
when k{rj^ — 77) <C 1 one has 

Q 1 

Avi = _e»kx-ifc»7, . ^ B^^ + hc 

4^3/2 ^5/2(^^_^)2 ^k + /i.C.. 

Hence, the resulting perturbations of the modulus have red power spectrum (1-4). 

The dependnce 6x1 oc (77* — ry)"^ is naturally interpreted in terms of the local 
shift of the "end time" parameter r/*. Indeed, with the background field given by 
(1-2), the sum Xc + ^Xi/V^i i-^-, the radial field including perturbations, is the 
linearized form of 

Xc[ry*(x) - ry] = — -i , (2-9) 

/i[?7*(x) -77] 

where 

77* (x) = 77* + (5r7*(x) (2-10) 

and 

So, the infrared radial modes modify the effective background by transforming the 
"end time" parameter 77* into random field that slowly varies in space, as given in 
eqs. (2-9), (2-10). 

It is worth noting that the infrared modes contribute both to the field (577* (x) 
itself and to its spatial derivative. The contribution of the modes which are super- 
horizon today, i.e., have momenta k < Hq, to the variance of the latter is given 

by 

{dir]4x)djr]^{x))k<Ho = ^ij ■ J^^^ y = Sij " "g^^ ^ ' (2-12) 

where A is the infrared cutoff which parametrizes our ignorance of the dynamics at 
the beginning of the conformal rolling stage. 
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§3. Effect of infrared radial modes on perturbations of phase: first 

order in h 

Let us sec how the interaction with the infrared radial modes affects the proper- 
ties of the perturbations of the phase 69. To this end, we consider perturbations of 
the imaginary part 6x2, whose wavelengths are much smaller than the scale of the 
spatial variation of the modulus (see ref.^^-* for details). Because of the separation 
of scales, perturbations 6x2 can still be treated in the linear approximation, but now 
in the background (2-9). 

Since our concern is the infrared part of 77* (x) , we make use of the spatial gradient 
expansion, consider, for the time being, a region near the origin and write 

77*(x) = 77*(0) - ViXi + ... , (3-1) 

where 

Vi = -5i?7*(x)|x=o , 

and dots denote higher order terms in x. Importantly, the field didjr}^{yi) has blue 
power spectrum, so the major effect of the infrared modes is accounted for by con- 
sidering the two terms written explicitly in (3-1). In this section we work at this, 
first order of the gradient expansion. Furthermore, we assume in what follows that 

|v| < 1 , (3-2) 

and in this section we neglect corrections of order v^. The expansion in |v| is 

legitimate, since the field v(x) has flat power spectrum, so the expansion in |v| is 
the expansion in /i, modulo infrared logarithms. We postpone to section 4.2 the 
analysis of the leading effect that occurs at the order v^. 

Keeping the two terms in (3T) only, we have, instead of eq. (1-2), 

= Mr?*(0)-r?-vx] ■ ^^'^^ 

This expression involves the combination 7?*(0) — (ry+vx). We interpret it as the local 
time shift and Lorentz boost of the original background (1-2). Note that the field 
(3-3) is a solution to the field equation (2-1) in our approximation. Our interpretation 
enables us to find the solutions to eq. (2-2) with the background (3-3) and initial 
condition (2-6): these are obtained by time translation and Lorentz boost of the 
original solution (2-4), (2-5): 

6xi~\K X, v) = e''i(''+^^)-^'?''*(0) • F{q, r?,(0) - r/ - vx) , (3-4) 

where the function F is still defined by (2-5), the Lorentz-boosted momentum is 

q = k-|-A;v, g = |q| = A; + kv , (3-5) 

and it is understood that terms of order must be neglected. We consider correc- 
tions of order 5^5^ 77* (x) and to this solution in sections 4.1 and 4.2, respectively. 
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We find from eqs. (3-3) and (3-4) that the perturbations of the phase again freeze 
out as A;[ry,,(x) — ry] — )■ 0, now at 

5e{^) = = ^ I ^-4^ ^^-^'^v*(^)A^ + h.c. , (3-6) 

This result implies that to the first order of the gradient expansion we limit ourselves 
in this section, the properties of the random field SO are the same as those of the 
linear field (2-7). Indeed, since ?7*(x) = 77*(0) — vx in (3-6), the infrared effects are 
removed by the field redefinition, 

iq = e-^*^^*W^ik, (3-7) 

where k and q are still related by (3-5). Due to the Lorentz-invariance of the measure 
d^k/k, the operators ^q, obey the standard commutational relations, while in 
our approximation, the field (3-6), written in terms of these operators, coincides with 
the linear field (2-7). We conclude that the infrared radial modes are, in fact, not 
particularly dangerous, as they do not grossly affect the properties of the field 69. 



§4. Statistical anisotropy 



4.1. First order in h 



To the first order in h, the non-trivial effect of the large wavelength perturbations 
Srj^: (x) on the perturbations of the phase occurs for the first time at the second order 
in the gradient expansion, i.e., at the order didjij^,. Let us concentrate on the effect 
of the modes of Sr]* whose present wavelengths exceed the present Hubble size. We 
are dealing with one realization of the random field (Jry*, hence at the second order 
of the gradient expansion, didjrj^ is merely a tensor, constant throughout the visible 
Universe. In this section we calculate the statistical anisotropy associated with this 
tensor. 

To this end, we make use of the perturbation theory in didjij^,. The background 
(2-9) is no longer the solution to the field equation (2-1) at the second order in the 
gradient expansion. The relevant combination entering eq. (2-2) for the perturbations 
of the imaginary part is now given by (see ref.^^^ for details) 

2h^x'^ = - + - . (4-1) 

The solution to eq. (2-2) with background (4-1) and initial condition (2-6) in the 
late-time regime is 

(-) ^gikx-ifcr;* (x) / ^ j^.j^. ^ \ 

The two non-trivial terms in parenthesis give the correction to the power spectrum 
of the phase perturbations due to the radial modes whose whavelengths exceed the 



Properties of scalar perturbations. 



9 



present Hubble size. The same correction is characteristic of the adiabatic pertur- 
bations, so we have finally 



(4-2) 



Notably, the angular average of the correction vanishes, so we are dealing with the 
statistical anisotropy proper. 

Neither the magnitude nor the exact form of the tensor didjTj^, — {l/3)Sijdkdkri.^ 
can be unambiguously predicted because of the cosmic variance. To estimate the 
strength of the statistical anisotropy, let us consider the variance 



where the notation reflects the fact that we take into account only those modes whose 
present wavelengths exceed the present Hubble size. In this way we arrive at the first 
non-trivial term in (1-5). Higher orders in the gradient expansion give contributions 
to the statistical anisotropy which are suppressed by extra factors of Ho/k. 

4.2. Order h?". contribution of deep infrared modes 

Let us now turn to the statistical anisotropy at the second order in h. Since 
the overall time shift r/*(0) is irrelevant, the major contribution at this order is 
proportional to the log-enhanced combination dirj^djij^ = ViVj. Hence, we use the 
two terms of the derivative expansion, explicitly written in (3T). 

To order v^, the function (3-3) is no longer a solution to the field equation (2-1). 
Instead, the solution is 

^7h*(0) vx] 

So, it is appropriate to study the solutions to eq. (2-2) in this background. It is 
straightforward to see that the solution that obeys the initial condition (2-6) is 

(5x^"^(k,x,r/) = e'9in'(^ll+''^)+^q^^^-'9T"?*(0) • F [q,-r{r)*{0) - r/ - vx)] , (4-4) 

where the Lorentz-boosted momentum is, as usual, g|| = 7(A;|| + kv), = k-^, 
q = ■^/{k + k\\v), J = (1 — v^)"^/^, and notations || and T refer to components parallel 
and normal to v, respectively. 

According to the scenario discussed in this paper, the phase perturbations 69 
freeze out at the hypersurface 77 = ?7*(0) — vx = r]^{0) — vx^^ and then stay constant in 
the cosmic time rj. It follows from (4-3) and (4-4) that at late times, the perturbations 
of the phase are 

5^(x) = / ^-4^ e*^"'«lFN+^'i^"^iq + h.c. , (4-5) 

where the operators are defined by (3-7) and obey the standard commutational 
relations, and we have omitted irrelevant constant phase factor. We see that to the 
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order v^, the effect of the deep infrared modes is encoded in the factor 7"^ in the 
first term in the exponent: the momentum of perturbation labeled by q is actually 
equal to p = (7^f/||, q^) • Accordingly, the power spectrum (omitting the correction 
discussed in section 4.1) is given by 

Hence, we have arrived at the last term in (1-5). Again, neither the direction of v 
nor its length can be unambiguously calculated because of cosmic variance; recall, 
however, that the value of |v|, and hence the parameter C2 in (1-5), is logarithmically 
enhanced due to the infrared effects, see (2-12). 

In the case of the curvaton mechanism of conversion of the phase perturbations 
into the adiabatic ones, both terms in the statistical ansotropy (1-5) are small because 
of the bound (2-8). On the other hand, the effect may be stronger in the case of the 
modulated decay mechanism. 



§5. Non-Gaussianity 

In this section we discuss non-Gaussianity of the phase perturbations (see ref.^^^ 
for details). Due to the symmetry 9 — >■ —9, the bispectrum {696969) as well as odd 
higher order correlators vanish identically and we have to deal with the trispectrum. 
Since the infrared modes are irrelevant to the leading order in h, the gradient ex- 
pansion of the modulus is useless, and we make use of the conventional perturbation 
theory. It is convenient to perform the calculation in terms of the perturbations of 
the modulus and phase. 



X 



6p\ / 69 



and employ the IN- IN formalism (see, e.g., ref.^^^). The late-time expectation value 
of an operator O is given by 



(O) 



Texp I i J dx^d^xHi 
■00 



{I) 



—00 



(5-1) 



where the time variable is x = 77 — 77* (note that x < 0) , T and T denote anti-time 
ordering and time ordering, respectively, and l-Lj and 0(7) are the interaction part 
of the Hamiltonian and operator O in the interaction picture. In what follows we 
are interested in the first nontrivial order in h, so the relevant part of T-Lj is 



ni = 



Xc6p 
V2 



{id,69f - {69T) 



(5-2) 



Upon substituting (5-2) into (5-1) one finds the following expression for the 
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connected part of the four-point function 



{de^dOySeje^) = j^j dx'QdxQd^x'd^x''xc{x'Q)xc{xQ) [Qx'^yQx"x^^QDp{x' ,x")- 

—oo 

- lm.Dp{x',x") {0{Xq - Xo)Q:r/xy^a;"zw " &{Xq - Xo)Ra;'ycyQx"z^v}] + 

+[y^z] + [y^w] (5-3) 



where 



Qx'xy = Im 



^-i^(x',x,x° = 0)^i5(x',y,2/° = 0) 



dx' 



^-Z)(x',x,x° = 0)-A_Dix',y,y' = 0) 



dx 



(5-4) 



D{x,y) and Dp{x,y) arc the pairings of (56' and 6p, respectively, D{x,y) = {66{x)d9{y)), 
etc. Recalling that in the interaction picture 69 and dp coincide with 6x2 /Xc and 
6xi, respectively, and using the results of sec. 2 one gets 



D{x',^,xo = 0) = ^^Jl I ^{-kx',f^H^%i-kx',)e^^(-'--^ 



(5-5) 



Dp{x,y) = J d^k^XQyQHfJ^{-kxQ)Hfl'^{-kyQy 



8k(x-y) 



(5-6) 



It is worth noting that in the derivation of eq. (5-3) we used the following property 
of the pairings 

D{x, y) = D*{y, x) , Dp{x, y) = D*p{y, x) . 

To find the explicit form of the functions Rx'xy and Qx'xy we substitute (5-5) 
into (5-4). Using the explicit form of the Hankel function, we get 

^ 8(27r)6 J kfkl «lfci,fc2,g, xo) 

8(27r)6 J k\kl ^ Q{ki,k2,q, x,) (5 7) 



x'xy 



where 



Riki,k2,q,0 = cos(C(fei + k2)) ■ [ekik2{q^ ~ {ki + fes)') + kf + kl - q^] + 
+ sin(C(fei + k2)) ■ aki + k2){kl + - q") 



Q{ki, k2,q,0 = sm(e(A;i + k2)) ■ [ekik2iq^ - (fci + fcs)') + kj + kj - q^] - 

- cos(e(A;i + k2)) ■ ^{ki + k2)ikl + k^ - q^), (5-8) 
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and 

q = ki + k2. 

The formula (5-3), converted to momentum representation, gives the general expres- 
sion for the trispectrum of the phase. 

Let us now consider the folded configuration, namely, the trispectrum in mo- 
mentum representation 

G(^) (Q, P, q) = (59 (I + q) 59 (| - q) <50 (-| + p) 59 (-^-^)) (5-9) 

with 

Qr^P:>q. (5-10) 
In real space, this limit corresponds to the following configuration 

|x — y| ~ |z — w| <C |x — z| ~ |y — w| , 

i.e., the four points x, y, z,w form, crudely speaking, an elongated parallelogram 
with the short edge e ~ |x — y| and the long edge L ~ |x — z|. 

To proceed, we consider in the small q limit the contribution to G*^^^ coming 
from the first term in the square brackets in (5-3). Making use of eqs. (5-6), (5-7) 
and (5-8), we find 

GW(Q,P,q)= ^l'',o,o\ o X (5-11) 

oo 

X J^^Q{k+,k-,q,OQ{p+,p-,q,x) [n/2 (90^5/2 (gx) + J5/2{qOJ5/2{qx)] , 



where 

k± = q/2±Q, p± = -q/2±Q, (5-12) 

and J5/2, ^5/2 are the Bessel functions of the first and the second kind, respectively. 
We see that the integral over ^ and % splits into a sum of products of two integrals 

00 

J{k+,k^,q)= J ^Qik+,k^,q,OJ5/2iqO (5-13) 



00 

yik+,k^,q) = J ^Q{k+,k-,q,OY,/2iqO ■ (5-14) 


These integrals converge in both the lower and upper limits*). This means, m par- 
ticular, that in the small q limit (if one considers k± and q as independent variables) 

To sec this ouc notes that the leading behaviour of Q{(,) at small ^ is ^'^j so the integrands 
are regular at 5 = 0. In the upper limit, the integrals become convergent if we make a rotation 
x'o x'Q+ie\x'o \ in eq. (5-3). This corresponds to the standard vacuum of the interacting theory. Note 
that we have to modify the contour of integration before splitting the expressions in the integrand 
into real and imaginary parts. Therefore, the expressions (5T3), (5T4) involve only decreasing 
exponents. 
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the leading behaviour of the integrals is 

Thus, we can neglet J. The direct evoluation of the integral (5-14) gives 
y{k^, k^,q) = -^yi^ [3(4 - klf - 2(4 + kl)q' - q'] = 




27r ^2 



7^ 



l-3(^) +0(gVQ'^ 



(5-15) 



where in the second line we make use of (5-12). We see that the contributions 
proportional to q~^ and q~^ vanish, as should be the case in view of the results of 
section 3. 

Upon substituting eq. (5-15) into eq. (5-11) we finally get in the small q limit 
the contribution to G*^^^ coming from the first term in the square brackets of (5-3): 



GW(Q,P,q) 



2^2 



h'TT 



1 



16 gg4p4 



1 



qP 

qP 



+ 0{q/Q). (5-16) 



This is actually the final expression for the trispectrum in the folded limit. To see 
this wc, first note that q is the momentum carried by the pairing Dp in (5-3), i.e. the 
regime (5-10) corresponds to the infrared effects associated with the perturbations 
of the modulus. It follows from (5-6) that at small momenta the leading behaviour 



of Dp is 



ReDp - 



ImDp - q^ 



Hence, only the first term in the square brackets in eq. (5-3) contributes to G^^^: 
the time integral of the second term as well as the second term itself are regular as 
g ^ 0. The crossing term (last line in (5-3)) does not contribute to (5-16) as well. 
The reason is that the momentum carried by Dp in that case is 2(Q + P) rather 
then q, so it is large. 



§6. Conclusion 



We conclude by making a few remarks. 

First, our mechanism of the generation of the adiabatic perturbations can work in 
any cosmological scenario that solves the horizon problem of the hot Big Bang theory, 
including inflation, bouncing/cyclic scenario, pre-Big Bang, etc. In some of these 
scenarios (e.g., bouncing Universe), the assumption that the phase perturbations 
are superhorizon in conventional sense by the end of the conformal rolling stage may 
be non-trivial. It would be of interest to study also the opposite case, in which the 
phase evolves for some time after the end of conformal rolling. 

Second, we concentrated in the first part of this paper on the effect of infrared 
radial modes, and employed the derivative expansion. The expressions like (3-6), 
which we obtained in this way, must be used with caution, however. Bold usage 
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of (3-6) would yield, e.g., non-vanishing equal-time commutator [0(x), 6'(y)], which 
would obviously be a wrong result. The point is that the formula (3-6) is valid in 
the approximation v = const; with this understanding, the equal-time commutator 
vanishes, as it should. 

Finally, the non-linearity of the field equation gives rise to the intrinsic non- 
Gaussianity of the phase perturbations and, as a result, adiabatic perturbations. The 
non-Gaussianity emerges at the order 0{h'^), and may be sizeable for large enough 
values of the coupling h. The form of the non-Gaussianity is rather peculiar in our 
scenario. Unlike in many other cases, the three-point correlation function vanishes, 
while the four-point correlation function of S6 (and hence of adiabatic perturbations) 
involves the two-point correlator of the independent Gaussian field (5??*. In view of 
the results of section 3, the correlation functions of 59 arc infrared-finite, at least to 
the order 0{h^). This is confirmed by the direct calculation in section 5. 

The authors are indebted to A. Barvinsky, S. Dubovsky, A. Frolov, D. Gor- 
bunov, E. Komatsu, V. Mukhanov, S. Mukohyama, M. Osipov, S. Ramazanov and 
A. Vikman for helpful discussions. We are grateful to the organizers of the Yukawa 
International Seminar "Gravity and Cosmology 2010" , where part of this work has 
been done, for hospitality. This work has been supported in part by Russian Foun- 
dation for Basic Research grant 08-02-00473, the Federal Agency for Sceince and 
Innovations under state contract 02.740.11.0244 and the grant of the President of 
Russian Federation NS-5525. 2010.2. The work of M.L. has been supported in part 
by Dynasty Foundation. 
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